Effects of Al doping on the structural and electronic properties of Mgi_a;Ala;B2 
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We have studied the structural and electronic properties of Mgj^_^ Ala;B2 within the Virtual 
Crystal Approximation (VGA) by means of first-principles total-energy calculations. Results for the 
lattice parameters, the electronic band structure, and the Fermi surface as a function of Al doping 
for < a; < 0.6 are presented. The ab initio VGA calculations are in excellent agreement with 
the experimentally observed change in the lattice parameters of Al doped MgBj. The calculations 
show that the Fermi surface associated with holes at the boron planes collapses gradually with 
aluminum doping and vanishes for x — 0.56. In addition, an abrupt topological change in the cr- 
band Fermi surface was found for x = 0.3. The calculated hole density correlates closely with existing 
experimental data for Tc{x), indicating that the observed loss of superconductivity in Mgj_^Ala;B2 
is a result of hole bands filling. 
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The discovery of superconductivity in the simple bi- 
nary lintermetallic compound MgBj with a as high as 
40 KEI has stimulated intense investigations, both from 
the experimental and the theoretical points of viewa. 
The superconducting transition temperature Tc for MgR* 
has been studied as a function of pressure and alloyingH. 
Pressure studies have showiLthat Tc decreases with ap- 
plied hydrostatic pressureffcl, which has been explained 
by an increaS|e|-pf the band filling of the boron cr-bands 
with pressurcQ'Ei. Thus, the change in the band filling in 
MgB2 under pressure is aa-effect of the reduction of the 
cell volume with pressureH'El. 

Experimentally it has been observed that the supercon- 
ducting transitkjM temperature of Mg]^_j.Alj;B2 decrease 
with Al dopingju, and superconductivity disappears for 
X > O.SQ'lI. According to band structure calculations 
of MgBj, electron doping-pduces the density of states 
(DOS) at the Fermi IcvcIEIIB. Based on the Rigid Band 
Approximation, An and PickettliJ analyzed the effect of 
Al doping on the DOS of MgB2, and found that the DOS 
at the Fermi level drops for x « 0.25. Structural char- 
acterization of Mg]^_2,Alj;B2 ^IaQws that the cell volume 
also decreases with Al dopingEll]. Therefore, the tr-band 
filling in Mg]^_j,Alj;B2 is expected to have two contribu- 
tions: the first is due to the electron doping, and the 
second is a result of the cell volume reduction as in the 
case of pressure effects. 

Measurements of the thermoelectric power, S, on 
Mg]^_^Ala;B2 for X < 0.1 have shown that the slope of 
the linear part of S{T) changes with Al doping, indicat,- 
ing changes in the Fermi surface due to electron dopingll3. 
More recently, it was shown tiiat the Raman spectra of 
Mg]^_^Al2;B2 for < X < O.WB show a pronounced fre- 
quency shift and a considerable change in the line- width 
for the i?2g phonon mode at a; w 0.3, which correlate,: 
a steeping in the behavior of Tc{x) with Al doping 
Although some of the observed effects of Al doping in 
Mg;^_^Alj,B2 can be interpretated qualitatively in terms 



of the Rigid Band Approximation as an effect of the 
cr-band filling, a quantitative analysis is essential in or- 
der to determine the interplay between electron doping 
and the structural, electronic, and transport properties 
of Mgi_,Al,B2. 

In this paper, we present a study of the effects of 
Al doping on the structural and electronic properties of 
Mg2_3,Ala;B2 for < a; < 0.6, using the ab initio Virtual 
Crystal Approximation. The calculated lattice para: 
eters are compared with available experimental dat. 
The evolution of the electronic band structure and the 
cr-band Fermi surface (FS) as a function of Al doping is 
analyzed. We correlate our results to the, experimentally 
observed behavior of Tc with Al dopingrlj. We show that 
the observed loss of superconductivity in Mg]^_j.Ala;B2 
can be explained by the filling of the hole bands. 

The Khon-Sham total energies were calculated self- 
consistently using the fuU-poteatial linearized augmented 
plane-wave mjethod (LAPW)ll3 as implemented in the 
WIEN97 codaiJ, where the core states are treated fully 
relativistically, and the semicore and valence states are 
computed in a scalar relativistic approximation. The 
exchange correlation potential was evaluated within the 
generalized gradient approximation (GGA), using the re- 
cent paranaeter-free GGA form by Perdew, Burke, and 
ErnzerholEj. We chose muffin-tin radii {Rmt) of 1.8 and 
1.5 a.u. for Mg and B, respectively, and used a plane- 
wave cutoff Rmt Km AX =9.0. Inside the atomic spheres 
the potential and charge density were expanded in crys- 
tal harmonics up to Z = 10. Gonvergence was assumed 
when the energy difference between the input and out- 
put charge densities was less than 1 x 10~^ Ryd. Special 
attention was paid to convergence of results by perform- 
ing the calculations for a sufficiently large number of k 
points in the irreducible wedge of the Brillouin zone for 
the omega structure (144 k points). The corrected tetp,- 
hedron method was used for Brillouin-zone integrationllj. 

The Al doping was modeled in the ah initio Virtual 
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Crystal Approximation (VCA)0S The Mg {Z = 12) 
sites are substituted by pseudo-atoms which have a frac- 
tional electronic charge {Z — 12 + x), depending on 
the Al concentration, x. This approximation is justified 
mainly by the fact that Al only has one electron more 
than Mg. The full-potential for the VGA system it is 
determined self-consistently foneach value of Al doping 
without shape approximationEj. The ab initio VGA as 
implemented in this work has been used very xecently 
to model G, Gu and Be substitutions in MgB2L3. The 
equilibrium lattice parameters were determined by total- 
energy calculations for each value of Al doping (a;=0.0, 
0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.5, and 0.6). Since the AIB2 
(omega) structure has two structural parameters (a and 
c) , we performed self-consistent total-energy calculations 
for nine different volumes and for nine different c/a ra- 
tios, in order to optimize both V and c/a for each Al 
concentration. For MgB2 we have obtained a = 3.083 
A and c = 3.526 A, which compares very well to tlie 
experimental values of a = 3.086 A and c ~ 3.524 AEI, 
respectively. 

In Fig. 1 we present the calculated lattice parame- 
ters (a and c) of Mg^_^Alj;B2 for < x < 0.6, and it 
can be seen that both a and & .decrease with Al dop- 
ing as observed experimentallyO'El. For comparison we 
have included the experimental data from Slusky et al.U, 
and we find that the change in a with Al doping is very 
well reproduced by the VGA calculations. Far the region 
0.1 < a; < 0.2 two values of c were reported^, which has 
been ascribed to the coexistence of two phases (Mg-rich 
and Al-rich phases). It is interesting to note that in this 
region the VGA values approximately reproduce the av- 
erage value. However, for x > 0.2 we find very good 
agreement between the experimental data and the VGA 
calculations (see Fig. 1). Although both cell parameters 
(a and c) decrease monotonically with Al doping, it is 
interesting to note that the slope for c as a function of Al 
doping is larger than for a. In order to understand this 
behavior in the context of bonding properties we have 
analyzed the change in the charge distribution with Al 
doping. Fig. 2(a) shows the charge density distribution 
in the (110) plane of MgB2. Mg nuclei are located at 
the corners of the map and B nuclei are at the (1/3, 1/2) 
and (2/3, 1/2) positions, all of them in the plane of the 
figure. We can see the directional, covalent B-B cr-bonds. 
In addition, there is a significant density of charge in the 
interstitial region giving rise to metallic-type bonding be- 
tween the Mg and B planes. The charge distribution and 
bonding properties of MgB2 have been calculated previ- 
ously and were discussed in detail in Ref. 19. Therefore, 
we concentrate on the influence of Al doping on the elec- 
tron density of MgB2. Fig. 2(b) shows the difference 
between the charge densities of Mgp 5AI0.5B2 and MgB2. 
It can be seen that charge transfer occurs from Al-B ions 
into the nearby interstitial region. We can see that the 
majority of this charge is distributed in the inter-plane 
region, and an important fraction of the charge is being 
transfered to the 7r-bond, while only a small fraction is at 



the cr-bond in the boron planes. This important increase 
in occupation of the 7r-bond with Al doping accounts for 
the strong decrease of the separation between planes (the 
c axis) of Mgj^_^Ala;B2. 

The evolution of the calculated DOS for Mg]^„2,Al2;B2, 
not shown here, shows a band broadening as a func- 
tion of Al doping, mainly as a consequence of the cell 
volume reduction. In addition, electron doping raises 
the Fermi level to higher energies. Both effects, band 
broadening and electron doping, contribute to reduce the 
density of states at the Fermi level in Al-doped MgB2. 
In Table I, we summarize the calculated lattice param- 
eters, cell volume, and the total density of states at the 
Fermi level [N{Ep)] for each of the studied Al concen- 
trations. We can see that N{Ef) decreases with Al dop- 
ing, from 0.72 for MgB2 to 0.26 states/eV per cell for 
MgQ4Alo.6B2. Therefore, in a BGS scenario this reduc- 
tion in N{Ef) accounts |for the decrease of Tc with Al 
doping in Mg^.^Al^Bjla. 

A careful analysis of the x-dependence of the electronic 
band structure, and in particular of the cr-band FS which 
has b cea s h own to be relevant for superconductivity in 
MgB2ByiH5m, provides a more detailed and quantita- 
tive description of the effects of Al doping. In Fig. 3 
we present the electronic band structure of MgB2. The 
(T-bands coming from the s-p boron orbitals, are strongly 
two-dimensional with very little dispersion along F-^, 
this dispersion can be characterized by the difference be- 
tween the Er and Ea energies (see Fig. 3). The Er and 
Ea energies correspond to the bottom and top of the tr- 
band in the T-A direction, respectively. The light-hole 
and heavy-hole cr-bands in MgBj form a FS consisting of 
two fluted cylinders surrounding the T-A line in the Bril- 
louin zone (see Fig. 3). The dependence of the energy of 
the cr-bands at F and A relative to Ep as a function of 
Al doping are shown in Fig. 4. We can see that both en- 
ergies, Er and Ea, decrease monotonically as a function 
of Al doping. More interestingly, the Fermi level reaches 
Er for X — 0.3 and Ea for x = 0.56. We find that the ra- 
dius of the cylinders decreases gradually with Al doping 
and at a critical concentration of a; = 0.3, the radius at 
kz — collapses and the FS takes the form of a sandglass. 
A three-dimensional view of the changes in the FS topol- 
ogy with Al doping are presented in Fig. 4. For x = 0.3 
the Fermi level in F is at a saddle point in the band struc- 
ture, and the transition through the saddle point results 
in the disruption of the neck, i.e_ the transition from a 
closed to an open section of FSEfl. For x > 0.3 the FS 
takes the form of two cones (see Fig. 4), and these finally 
vanishes at the second critical concentration {x ~ 0.56) 
when the hole bands have been filled. These changes in 
the hole FS are expected to be accompanied by various 
kinds of electronically driven anomalies, including lattice 
dynamics and transport properties. 

As was mentiaued above, Raman spectroscopy studies 
on Mg]^_^Ala;B2l3'Ll show a pronounced shift and a consid- 
erable change in the line- width of the i?2g phonon mode 
at a; « 0.3. Additionally, a steeping of the Tc decrease has 
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beea-,abserved at an Al concentration of approximately 
O.Soptl. These changes in both the structural and the 
superconducting properties seems to be strongly related 
to the abrupt change in the FS topology, which occurs 
for X = 0.3 (see Fig. 4). In order to establish a more 
direct comparison between the FS evolution and the su- 
perconducting properties with increasing the Al doping, 
we have calculated the hole FS area as a function of x, 
which is proportional to the hole density at the Fermi 
level. In Fig. 4(b) we compare the calculated normalized 
FS area, Ap s{x) / Ap s(Q) , with the normalized supercon- 
ducting critical temperaure, Tc(x) /TdO). The experi- 
mental data for Tc{x) were taken from Ref. 8. We can 
see that for the low concentration region {x < 0.25), be- 
fore the E2g phonon frequency shiftEi, the drop of Tc is 
directly related to the change in the hole density. This 
view is in agreement with recent resultsj-af NMR exper- 
iments on Al-doped MgB2 for x < O.lcJ. In the high 
concentration region {x > 0.25), the behavior of is 
determined by the FS area but.|lhc importance of the 
phonon-renormalization is clearcfcl. In this way, the FS 
area and follow the same behavior with Al doping in 
the whole range (0 < a; < 0.6), indicating a close relation 
between the changes in the cr-band FS and the loss of 
superconductivity in Al doped MgBj. 

In summary, we have performed a first-principles study 
of the effects of Al doping on the structural param- 
eters, the electronic structure, and the cr-band FS of 
Mg^_^Ala;B2, using the Virtual Crystal Approximation, 
(i) We find that the ab initio VGA calculations are in 
excellent agreement with the experimentally observed 
changes in the lattice parameters as a function of Al 
doping, (ii) The analysis of the charge density shows 
that an important portion of the Al-electrons are at the 
inter-plane region and only a small fraction at the B-B 
planes, providing an explanation for the strong change 
of the c-axis and the small change in the a-axis with the 
Al concentration, (iii) The hole FS gradually collapses 
with Al doping and vanishes for x = 0.56. An abrupt 
topological change was found for x = 0.3, which corre- 
lates with the frequency shift of the i?2g phonon mode 
and the steeping in the Tc{x) decrease. Additionally, the 
critical concentration of a; = 0.56 at which the hole FS 
disappears, corresponds to the experimentally observed 
Al concentration (0.5-0.6) for which Tc{x) vanishes, (iv) 
We find that the behavior of the calculated cj-band FS 
area with Al doping correlates with the superconduct- 
ing critical temperature Tc{x). Consequently, the ob- 
served loss of superconductivity in Mg]^_j.Ala;B2, can be 
explained as a result of the filling of the hole bands. 
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FIG. 1. Lattice parameters a and c for Mgj_^Al2,B2. Ex- 
perimental data from Ref. [6] (solid circle) and calculated 
using VGA (open circle). 



FIG. 2. Electronic charge density for the MgB2(110) plane 
(top) and the charge density difference between Mg^, 5AI0.5B2 
and MgBj (bottom). 



FIG. 3. Electronic band structure and hole Fermi surface 
for MgBj at the calculated lattice constants (see Table I). 



FIG. 4. (a)Energy position of the a-band at F (Er) and 
A (Ea) relative to Ef for Mgj_^Al2,B2. In the inset, the 
hole Fermi surface for x — 0.25, 0.3, and 0.35. (b) Calculated 
normalized holes FS area, A{x)/A{0) (solid line) and the nor- 
malized experimental values (Ref. [8]) of the superconducting 
critial temperature, Tc{x) /Tc{0) (open circles). 



TABLE I. Calculated lattice parameters, cell volume, and 
density of states at the Fermi level, N{Ef) in states/eV unit 
cell, for Mgj^_^Ali:B2 as a function of Al doping (x). 
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a(A) 


c(A) 




N{Ef) 


0.0 


3.083 


3.526 


29.02 


0.72 


0.1 


3.076 


3.486 


28.56 


0.68 


0.2 


3.072 


3.448 


28.18 


0.64 


0.25 


3.070 


3.424 


27.95 


0.60 


0.3 


3.063 


3.403 


27.65 


0.55 


0.35 


3.059 


3.386 


27.44 


0.48 


0.4 


3.055 


3.367 


27.21 


0.43 


0.5 


3.047 


3.338 


26.84 


0.33 


0.6 


3.039 


3.315 


26.51 


0.26 
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Figure 1 . O. de la Pena, et al. 




Figure 2. O. de la Pena et al. 
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